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, ^ ■ Abstract 

\^ • We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared 

^Sj \ flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm 

observations. Our goal is to understand the physics determining the infrared spectral energy distribution 
(SED). The ratio of the total infrared luminosity Xtir, to the star-formation rate (SFR) is tightly correlated 
by a power-law to specific SFR (SSFR), and Ltir is a good SFR indicator only for galaxies with the 
largest SSFR. Wc discovered a tight linear correlation for normal galaxies between the radiation field 
strength of dust heating, estimated by infrared SED fits [Uh)-, and that of galactic-scale infrared emission 
^ ^ (t^TiR oc Ltik/R^), where R is the optical size of a galaxy. The dispersion of Uh along this relation is 

I I 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the Uh/UTm ratio 

2 ■ can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust 

» I screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the 

. dust column density, supporting this interpretation. In the large C/tir limit, the data of circumnuclear 

I— I' starbursts indicate the existence of an upper limit on Uh, corresponding to the maximum SFR per gas 

mass of 10 Gyr~^. We find that the number of galaxies sharply drops when they become optically 
' thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by 

, the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. 

' Implications are discussed for the A'/ni-size relation, the Kennicutt-Schmidt relation, and galaxy formation 

I in the cosmological context. 

lO ' Key words: galaxies: formation — galaxies: ISM — galaxies: statistics — infrared: galaxies — ISM: dust, 

' extinction 

o 

■ 1. Introduction physical quantities (e.g.. Dale et al. 2001; Dale & Hclou 

; 2002; Draine et al. 2007; da Cunha et al. 2008, 2010; 

.5^ ■ Dust grains in interstellar medium (ISM) absorb ultra- Muhoz-Mateos et al. 2009; Malek et al. 2010). 

' violet and optical photons emitted from stars, and the The intensity of the cosmic infrared background radia- 
, energy is re-emitted from dust grains in infrared and sub- tion indicates that about half of the energy originally ra- 
_ C^_i millimeter wavelengths. The spectral energy distribution dialed from stars is eventually re-emitted from dust grains 
(SED) of this dust emission has a broad peak of the mod- (Hauser & Dwek 2001). Heavily obscured star formation 
ified blackbody in far-infrared bands from large grains activity is difficult to detect by observations in ultravi- 
in thermal equilibrium with temperatures of ~ 10-50 olet/optical observations, and infrared/submillimeter ob- 
K. The SED in near- to mid-infrared bands is, on the servations should crucially be important to get the full 
other hand, characterized by emissions from small grains picture of the formation and evolution of galaxies. In the 
and polycyclic aromatic hydrocarbons (PAHs) that have last decade, our observational knowledge about the hid- 
a wide distribution of temperature because of stochas- den side of galaxy evolution has significantly improved 
tic single-photon heating (Desert et al. 1990; see Draine thanks to the advanced infrared/submillimeter facilities 
2003 for a review). Physical dust models about composi- such as JCMT/SCUBA, Spitzer, AKARI, and Herschel. 
tions and size distribution have been developed to calcu- Moreover, the Atacama Large Millimeter/submillimeter 
late SEDs of dust emission for a given heating radiation Array (ALMA) will revolutionize our understanding of 
field strength, U, and one can predict the global SED in- galaxy evolution in the near future. In order to extract 
tegrated over a galaxy by summing up the emission from useful information from these data, we need good theoret- 
dust with various values of U. Such models have been ical understanding and models of galactic dust emission 
applied to fit observed SEDs of nearby galaxies to derive to be compared. 
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However, it is not straightforward to predict global (i.e., 
galactic-scale) SEDs of dust emission for a theoretical 
model galaxy, in contrast to direct stellar emission that 
can be calculated by the method of stellar population syn- 
thesis. This is because the physics determining the global 
SED of dust emission and its relations to other physical 
properties have remained highly elusive (see Walcher et 
al. 2011 for a review). An approach often taken in the 
literature (e.g., Guiderdoni et al. 1998; Takeuchi et al. 
2001; Chary & Elbaz 2001; Lagache, Dole, & Puget 2003; 
Valiante et al. 2009) is to relate the SED parameter (e.g., 
the dust temperature T^) to the total infrared luminosity, 
^TiR (defined as the bolometric luminosity of dust emis- 
sion), based on the observed correlation between the two 
(higher Tj, for larger Ltir; Soifer et al. 1987a; Soifer & 
Neugebauer 1991; Chapin et al. 2009). However, there 
is a large scatter along the mean L^ir-Tj, relation (see, 
e.g., Hwang et al. 2010 for recent data), and it is phys- 
ically unreasonable to relate an extensive quantity Ltir 
that scales with the system size to an intensive quantity 
Td that docs not (Totani & Takeuchi 2002). Extension of 
this relation including another parameter has also been 
discussed, e.g., galaxy size (Devereux 1987; PhiUips & 
Disney 1988; Lehnert & Heckman 1996; Chanial et al. 
2007; Rujopakarn et al. 2011a, b), but a consistent physi- 
cal picture has not yet been established. 

To make a completely theoretical prediction, one must 
solve transfer of dust heating radiation taking into ac- 
count complicated geometry in a galaxy. Infrared SED 
models by such ab initio approaches have also been devel- 
oped (Silva et al. 1998; Takagi et al. 2003a, b; Dopita et 
al. 2005), and applied to predict infrared SEDs in galaxy 
evolution models (Granato et al. 2000; Baugh et al. 2005; 
Fontanot et al. 2007). However, it is difficult to calculate 
detailed geometry within galaxies in models on the cos- 
mological scale, and one must add many uncertain and 
adjustable model parameters, making it difficult to ex- 
tract useful information from comparison with observed 
data. 

Therefore it is still important to investigate the relation 
between infrared SEDs and various physical properties to 
find any correlation that would be useful to better under- 
stand the galactic-scale radiation from dust grains. The 
recently released catalogs of the AKARI All Sky Survey 
(Ishihara et al. 2010; Yamamura et al. 2010) for sources 
detected by the AKARI satellite (Murakami et al. 2007) 
provide us with a new opportunity to investigate nearby 
galaxies for this purpose. The Infrared Camera (IRC, 
Onaka et al. 2007) and the Far-Infrared Surveyer (FIS, 
Kawada et al. 2007) on board the satellite have detected 
more than 800000 and 400000 sources in 9-18 and 65-160 
/im bands, respectively. For recent studies on galaxies us- 
ing samples based on these catalogs, see Takeuchi et al. 
(2010), Buat et al. (2011), Goto et al. (2011a, b), and 
Yuan et al. (2011). 

Especially, the four FIS bands at 65-160 ^m should be 
very useful to accurately measure the modified blackbody 
peak of the dust emission SEDs, as they cover the wave- 
length range around the peak by a larger number of pho- 



tometric bands with better flux sensitivity and angular 
resolution than those of the previous infrared all-sky sur- 
vey by the Infrared Astronomical Satellite (IRAS, Soifer 
ct al. 1987b). In this paper we perform SED fittings of 
physical dust models to the AKARI data of more than 
1600 low-redshift galaxies, and compare the SED param- 
eters with other physical properties of galaxies obtained 
by optical and 21-cm observations. The goal of this paper 
is to find the key physical quantities and laws to determine 
the global SEDs of dust emission (especially around the 
thermal peak) , and give physical interpretations to them. 

In §2 we describe the construction of the galaxy sam- 
ples used in this work. The physical dust models and SED 
fittings to the AKARI data are described in §3. We then 
present theoretical background and considerations in §4, 
which will be used to interpret the results obtained in this 
work. The main results are presented in §5, followed by 
discussions in §6 and summary in §7. We adopt the stan- 
dard cosmological parameters of Hq = 70 kms~^ Mpc~^, 
= 0.3, and f^A = 0.7, though the choice of the cos- 
mological parameters hardly affects our results based on 
galaxies at low redshifts. 

2. The Samples 

2.1. The AKARI All-Sky Survey 

We use the AKARI/IRC All-Sky Survey Point Source 
Catalog (PSC, ver. 1.0) and the AKARI/FIS All-Sky 
Survey Bright Source Catalog (BSC, ver. 1.0), which 
were released in March 2010 including 870973 and 427071 
sources, respectively. The IRC provides fluxes in two 
photometric bands, SOW and L18W, whose central wave- 
lengths (Ac) are 9 and 18 /im, respectively. The FIS pro- 
vides fluxes in the four photometric bands of N60, WIDE- 
S, WIDE-L, and N160, whose central wavelengths are 65, 
90, 140, and 160 ^m, respectively. The sensitivity limits 
corresponding to 50% detection completeness are shown 
in Table 1, and the relative spectral response functions 
(RSRFs) for the six bands are shown in Fig. 1. See the 
release notes of these catalogs available on the AKARI 
web site for more details. 

Since the SED fltting is the central issue of this work, 
treatment of flux errors is important. In this work we con- 
sider only statistical flux errors associated with each flux 
measurements, and systematic errors that are common 
for all sources (e.g., flux calibration errors) are not taken 
into account. Systematic errors would result in systematic 
bias in SED fit parameters (e.g., a systematic shift of dust 
temperature) , but it does not seriously affect the study of 
correlation or trend between the flt parameters and other 
galaxy properties. Statistical errors can be divided into 
two parts: one is the error that does not depend on source 
fluxes (denoted as AFconst), and the other is the error that 
is proportional to flux {AF = eF) . The total flux error is 
assigned for a flux measurement in a band as: 

AFtot= [(A^^eonst)' + (e^^)']^ . (1) 

For faint sources Ai^const is dominant while eF is domi- 
nant for bright sources. The la values of AFconst and e 
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Fig. 1. The relative spectral response function (RSRF) 
against energy flux for the two AKARI/IRC bands {S9W 
and L18W) and the four AKARI/FIS bands {N60, WIDE-S, 
WIDE-L, and N160). The RSRFs are normalized by their 
maximum values. 

for each band are estimated from the signal-to-noise dis- 
tributions presented in the release notes of the AKARI 
catalogs, and they are summarized in Table 1. 

2.2. Cross-Matches with Optical and 21cm Data 

We construct the following three samples by cross- 
matching the AKARI catalogs with the existing catalogs 
of nearby galaxies selected in other wavelengths. The first 
sample is the cross-match with the galaxies detected by 
the Sloan Degital Sky Survey (SDSS), and the second with 
the Hi Parkes All-Sky Survey (HIPASS) selected by 21- 
cm flux. We regard these as the two main samples in 
this paper, and we can check the dependence on the sam- 
ple selection by these two. In addition, we construct the 
third sample to extend the dynamic range of surface den- 
sity of star formation rate (SFR), which is the cross- match 
with the Kennicutt (1998, hereafter K98) sample of nearby 
circumnuclear starburst (SB) galaxies. The summary of 
constructing these samples, showing galaxy numbers at 
various stages of selections, is presented in Table 2. 

2.2.1. The AKARI-SDSS Sample 

We first extracted a sample of 1 145 617 SDSS galax- 
ies having measured redshifts from the DR7 NYU Value- 
Added Galaxy Catalog (NYU-VAGC, Blanton et al. 
2005). The information of apparent Petrosian magnitudes 
in the five ugriz bands, (Galactic) extinction- and K- 
corrcctcd absolute magnitudes, and r-band Petrosian 90% 
angular radii that include 90% of Petrosian fiux {Opqq) will 
be used in the analyses of this work. We take 9pqo from 
several size estimations available for SDSS galaxies, since 
this size is not sensitive to central concentration (Nair et 
al. 2010). The angular sizes are converted into physi- 
cal radii as Rpgo = dA{z) ^pgo, where dA is the standard 
angular diameter distance to redshift z. 

These galaxies were cross-matched separately with the 
AKARI/IRC and AKARI/FIS catalogs, with the match- 
ing radius of 5 and 20 arcsec, respectively. The mean 
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positional errors of sources in these catalogs are 0.8 and 
6 arcsec, and essentially all objects should fall within the 
matching radii. [See also Table 1 for the instrumental 
point spread functions (PSF) in each band.] To avoid 
contamination of random associations, we removed any 
duplications, i.e., AKARI galaxies having more than two 
SDSS counterparts, and vice versa. We then obtained 390 
and 6500 matches of SDSS galaxies with the IRC and FIS 
catalogs, respectively. The expected random associations 
arc 25.1 and 143.7 for SDSS-IRC and SDSS-FIS, if the 
sources arc uniformly distributed on the sky. These num- 
bers are sufficiently small for studying the overall and sta- 
tistical trends of galactic properties. Because the number 
of cross-matches with the IRC catalog is small, we con- 
struct the AKARI-SDSS sample based on the correlation 
between the SDSS and FIS catalogs, without requiring 
detection by IRC. 

To keep good quality of the data, we further selected 
galaxies satisfying the following conditions: z > 0.005, r < 
16.0 mag, 6*^90 > 6.0", and the AKARI/FIS flux quality 
flagi FQUAL > 3 in at least two FIS bands. The redshift 
condition is to remove too close objects whose redshift- 
based distance is rather uncertain by the peculiar velocity 
effect. The magnitude and size conditions are to keep 
good quality of size estimates and to make the seeing effect 
negligible. The condition for FQUAL is for reliable SED 
flttings. Then 878 galaxies qualified as the final AKARI- 
SDSS sample. The numbers of galaxies with FQUAL > 1 
or 3 in each FIS band are shown in Table 3. This table 
shows that the 90- and 140-^m bands are more sensitive 
than the other two bands and most galaxies are detected 
in these bands with FQUAL > 3. 

To these galaxies we added physical quantities derived 
from SDSS spectra by the MPA/ JHU group when they are 
available, by cross-matching the spectroscopic ID [plate 
number, Mean Julian Date (MJD), and fiber number]. 
Stellar mass (M*, Kauffmann et al. 2003; Salim et al. 
2007), SFR (V;, Brinchmann et al. 2004), and gas-phase 
metallicity (Tremonti et al. 2004) of the DR7 version are 
available for 550, 567, and 246 galaxies, respectively, in 
the total sample of the 878 galaxies. The metallicity is 
given in Aq = 12 -f logio(0/H) in the original data, but 
in this work we use metallicity Z converted^ from 
using AoQ = 8.69 and Zq = 0.02. We also added the V 
band optical depth Ty of stellar light attenuation by dust 
in galaxies, estimated by Tojeiro et al. (2009) using the 
VESPA algorithm (the one-parameter dust model in the 
RunID=l catalog), to 534 galaxies. 

Wc also added information of neutral hydrogen (Hi) 
mass from 21-cm observations in the literature. SDSS 
galaxies were cross-matched with the optical counterpart 

^ See the release note of the AKARI/FIS catalog for the definition 
of FQUAL. It is recommended not to use the flux data when 
FQUAL < 2 for a reliable scientiflc analysis. 

^ There are systematic uncertainties in estimates of Aq, Z, and 
their solar values (Grevesse et al. 2010; Moustakas et al. 2010). 
We will mostly discuss correlations between metallicity and 
other physical parameters, and absolute uncertainty does not 
signiflcantly affect our conclusions. 
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Table 1. Summary of the AKARI All-Sky Survey Parameters 



instrument 




IRC 






FIS 




band name 


S9W 


L18W 


N60 


WIDE-S 


WIDE-L 


N160 


central wavelength Ac [/^m] 


9 


18 


65 


90 


140 


160 


flux sensitivity* [Jy] 


0.09 


0.17 


2.6 


0.39 


2.9 


5.9 


AFconst [Jy] 


0.02 


0.04 


0.48 


0.11 


0.28 


1.26 


. [%] 


5.6 


6.7 


10 


10 


10 


10 


PSF [arcsec FWHM] 


5.6 


5.6 


37 


39 


58 


61 



*The flux where the source detection completeness is 50%. 



Table 2. The Summary of Sample Construction 



AKARI-SDSS 


AKARI-HIPASS AKARI-SB 


original catalog 


1145617 


4315 


36 


match with AKARI/FIS 


6500 


1527 


34 


FIS flux quality flag 


1138 


840 


33 


distance cut* 


996 


711 


32 


magnitude & size cuft 


878 




31 


.f'TIR,K98/^TIR > 0.5 






24 


final sample 


878 


711 


24 


stellar massj 


550 








star formation ratef 


567 








metallicityl 


246 








TV 


534 








H I gas mass 


257 


711 





H2 gas mass 








24 



Notes. — The upper six rows show the numbers of galaxies after applying certain selection criteria in the 
process of constructing the samples. The lower six rows show the numbers of galaxies for which a certain 
physical quantity is available. 

*The condition oi z> 0.005 for AKARI-SDSS and AKARI-HIPASS, while z > 0.005 or availability of rcdshift- 
independent distance is required for AKARI-SB. 

^The condition of r < 16.0 mag and OpgQ > 6" is adopted for AKARI-SDSS, while a galaxy without diameter 
size in the K98 sample is excluded from AKARI-SB. 
^Physical quantities derived from SDSS optical spectra. 



positions of 21-cni-detected galaxies in the data sets of 
Springob et al. (2005), the ALFALFA survey (Giovanehi 
et al. 2007; Kent et al. 2008; Stierwalt et al. 2009), the 
Equatorial Survey (Garcia- Appadoo et al. 2009; West 
et al. 2010), and the GALEX Arecibo SDSS Survey 
(Catinella et al. 2010). We found 257 Hi flux data avail- 
able for the AKARI-SDSS sample, and calculated the H I 
mass from the standard formula'^: 



2.36 X 10^ 



+ \ Mpc / VJykm/s 



S\, 



(2) 



where Sint = J Si,dV is the velocity- integrated flux den- 
sity and c?i is the standard luminosity distance. For com- 
parison, we also constructed the SDSS-Hl sample, for all 
SDSS galaxies with available H I data without requiring 
the detection by AKARI. Other selection criteria were 
kept same as the AKARLSDSS sample. The SDSS-Hl 
sample consists of 2413 galaxies. 



^ The factor of (I-I-2) ^ arises as follows. The Hi mass can be 
written as Mm (X J S^du = d^ J (du/dV) Sy dV, where 1^ is 
the observed frequency. In Hi observations, the velocity V is 
defined as u = !^2i/(l -I- V/c) where 1/21 = 1.42 GHz, and hence 
du/dV oc (1 + V/c)~^ ~ (1 -I- when the intrinsic velocity 

dispersion AV in a galaxy is much smaller than cz. 



Table 3. The Number of Galaxies against AKARI/FIS Flux 
Quality Flags 









FIS bands [fJ-m] 




sample 


condition 


65 


90 


140 


160 


AKARI-SDSS 


EQUAL > 1 


878 


878 


878 


878 




EQUAL > 3 


178 


876 


851 


136 


AKARI-HIPASS 


EQUAL > 1 


710 


711 


711 


710 




EQUAL > 3 


173 


704 


700 


239 


AKARI-SB 


EQUAL > 1 


24 


24 


24 


24 




EQUAL > 3 


21 


22 


24 


23 



2.2.2. The AKARI-HIPASS Sample 

Since the gas amount of galaxies is crucial for the dust 
amount, we also cross-correlated the AKARI galaxies with 
the HIPASS catalog (HICAT, Meyer et al. 2004), which 
is a catalog of 4315 galaxies selected by the Hi 21-cm line 
flux covering the entire southern sky. We selected HIPASS 
galaxies having optical counterparts from the HOPCAT 
catalog (Doyle et al. 2005), for which the information of 
21-cm line flux, 21-cm redshift, optical magnitudes of Bj, 
R, and /, and the major axis size Oa.hop of i?-band im- 
ages are extracted from the catalog. The distances were 
calculated from 21-cm redshifts, and then Mhi was calcu- 
lated from eq. (2). Absolute optical magnitudes were also 
calculated, correcting the Galactic extinction but ignoring 
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Fig. 2. Correlation between the r-band 90% Petrosian an- 
gular radius (9pgo) of SDSS and the il-band major axis size 
(^A,HOp) of HOPCAT for galaxies commonly included in the 
SDSS and HIPASS samples. Filled circles arc galaxies de- 
tected by AKARI/FIS. 

ii'-correction. 

Then the HIPASS galaxies were cross-matched with the 
AKARI catalogs by the same procedures as those for the 
AKARI-SDSS sample. We found 161 and 1527 matches 
with the IRC and FIS catalogs, respectively. Following the 
AKARI-SDSS sample, we do not require detection by IRC 
for the AKARI-HIPASS sample. To keep good qualities, 
we further require the same criteria about redshift [z > 
0.005) and FIS quality flags {EQUAL > 3 at least in two 
FIS bands) as the AKARI-SDSS sample, resulting in the 
final AKARI-HIPASS sample of 711 galaxies. 

Because the optical magnitude and size of HOPCAT 
are from photographic plates, their quality and reliabil- 
ity should carefully be examined. To directly compare 
the AKARI-SDSS and AKARTHIPASS samples, we con- 
verted the HOPCAT magnitude (i?) and size (6'a.hop) 
into those in the SDSS catalog (r magnitude and ^pgo), 
and will use these in the following analyses. For this con- 
version, we first searched for galaxies that are common 
in SDSS and HOPCAT by optical positions, and found 
103 galaxies that satisfy the same criteria as the AKARI- 
SDSS galaxies about redshift, magnitude, and size. In 
these galaxies 47 are detected by AKARI/FIS. From com- 
parison between SDSS r and HOPCAT R bands, we apply 
the following conversion: i?(HOP) - r(SDSS) = 0.09, and 
the dispersion around this relation is 0.46 mag. Figure 
2 shows the comparison between the SDSS size Opgo and 
the HOPCAT size 0a,hop, and the ratio 0^,hop/^P9O sys- 
tematically varies with 0a. hop- We therefore convert the 
HOPCAT size into Opqq by the following equation: 

/ ^A.HOP 



log 



10 



-0.489 + 0.245 logio f^^S ) (3) 



\ arcsec 



which is shown in Fig. 2 by the solid line. 




/ 



Fig. 3. The histogram of Ltir,K98/^tir for the K98 star- 
burst sample, where Ltir,K98 is the total infrared lumi- 
nosity from the circumnuclear starburst regions given in 
Kennicutt (1998) and -Ltir is the total infrared luminos- 
ity from the whole region of galaxies estimated by SED fit- 
tings to the AKARI data in this work. Only galaxies having 
^TiR.KQs/^TiR > 0.5 are used as the AKARI-SB sample. 



2.2.3. The AKARI-SB Sample 

The AKARI/FIS catalog was cross- matched with the 36 
galaxies of the K98 circumnuclear starburst sample by the 
same procedures as the two samples described above, and 
we found that 34 galaxies are detected by FIS. To keep 
the quality of the data, we again adopt the same condi- 
tion about the FIS flux quality flags. For these galaxies 
we collected the redshift-independcnt distance measure- 
ments from the NASA Extragalactic Data base (NED), 
and adopted the mean value when more than one mea- 
surements are available for a galaxy. For galaxies that do 
not have a redshift-independent distance measurement, we 
calculated the distances from redshifts, but NGC 5194 was 
removed from the sample because z < 0.005. We will use 
the information of the angular diameter size 9o , SFR sur- 
face density Ssfr, and H2 gas mass surface density Sh2 
given in K98. Here, NGC 3256 was removed because 9d 
is unavailable. 

In K98, SFR and H2 mass were estimated from Ltir 
and CO line luminosity, respectively. As argued in K98, 
Ltir is a good SFR indicator for these starburst galax- 
ies because Ltir is dominated by young stars and dust 
opacity is large. The SFR-Ltir conversion factor of K98 
assumes the Salpeter (1955) initial mass function (IMF), 
and we divided K98 SFRs by 1.5 to match the SFRs in the 
AKARI-SDSS sample assuming the Kroupa (2001) IMF 
(Brinchmann et al. 2004). In this case the SFR (denoted 
as ?/;) is related to Ltir as: 



Ltir 



8.7 X 10^ 



(4) 



Lq \Mq yr' 

It should be noted that in this work this relation is used 
only for the AKARI-SB sample; we do not use Ltir as a 
SFR indicator for the AKARI-SDSS or AKARI-HIPASS 
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Fig. 6. Comparison between the AKARI-SDSS and 
AKARI-HIPASS samples in the plot of Pctrosian 90% 
angular radius (6p9o) versus SDSS r magnitude, after the 
conversion of the HOPCAT R and Sa.hop into the SDSS 
quantities as described in the text. 

samples, because a significant fraction of L^ir comes from 
dust heating by relatively aged stars (see §5.1). SFRs from 
SDSS spectra will be used for the AKARI-SDSS sample. 

The surface densities (Sh2 and Ssfr) given in K98 
are converted into the total quantities (Mh2 and ip)i o.g., 
Mh2 = ttR^ Sh2, where i? = 6*0 o?a/2 is the physical radius 
of circumnuclear starburst regions. Note that the "total" 
quantities in the AKARI-SB sample are the total within 
the circumnuclear starburst regions rather than the entire 
galaxies. Since the AKARI photometries are for the whole 
region of galaxies, we must examine whether the AKARI 
fluxes are mostly coming from the circumnuclear regions. 
For this purpose, we compared the total infrared lumi- 
nosity from a circumnuclear region, -/jtir,K98j of the K98 
sample (calculated from Esfr and 9d) with Ltir for the 
entire galaxy newly obtained by our SED fit to the AKARI 
data (see §3), in Fig. 3. It can be seen that many galaxies 
have Ltir ixiR.KQS) indicating that the total infrared 
luminosities of galaxies are dominated by those from cir- 
cumnuclear starburst regions. We adopted the condition 
of £tir.K98 1^ 0.5 -Ltir for our sample, to remove galaxies 
whose infrared luminosity is mostly emitted from outer 
regions. Then finally we obtained the AKARI-SB sample 
consisting of 24 galaxies. 

2.3. The Basic Properties of the Sample Galaxies 

Figure 4 shows redshift versus the AKARI 90-/im flux 
for the final AKARI-SDSS and AKARI-HIPASS samples. 
The sensitivity of the 90-/zm band is the best among the 
four FIS bands. Sources are distributed down to the flux 
sensitivity limit, but the mean flux is considerably brighter 
than the flux distribution of all FIS-BSC sources reported 
in the release note, because we have set the condition of 
FQUAL > 3 at least in two FIS bands. The mean redshift 



of the AKARI-SDSS sample {z - 0.03) is slightly higher 
than that of the AKARI-HIPASS sample (z - 0.01). 

Figure 5 shows angular size versus apparent optical 
magnitude (r for SDSS and R for HIPASS). Compared 
with the typical SDSS magnitudes, the AKARI-SDSS 
galaxies are much brighter, indicating that the AKARI 
sensitivity limit is brighter than SDSS for typical galaxy 
SEDs. In Fig. 6, we show the size-magnitude plot of 
the AKARI-HIPASS galaxies when the R magnitude and 
Oa.hop are converted into r and 0_p9o, in comparison with 
the AKARI-SDSS galaxies. The distributions of the two 
samples are similar, indicating that the conversion of mag- 
nitude and size is working well. 

Figure 7 shows the optical color versus absolute magni- 
tude relation (see also Takeuchi et al. 2010). Our samples 
cover a wide range of galaxy colors, though the AKARI- 
SDSS galaxies have a slightly redder mean color than 
the mean SDSS distribution, indicating that we are sam- 
pling dusty galaxies. Red and brightest SDSS galaxies are 
not detected by AKARI, because they are early-type red- 
sequence galaxies with no or low star formation activity. 

Figure 8 shows the stellar mass (M*) versus metallic- 
ity relation of the AKARI-SDSS sample, which indicates 
that we are sampling relatively massive and high metallic- 
ity galaxies compared with general SDSS galaxies. This is 
reasonable for our samples limited by the AKARI sensi- 
tivity, because massive star-forming galaxies are generally 
dusty and high metallicity is required for a large amount 
of dust. 

Figure 9 shows V band attenuation Ty versus specific 
star formation rate (SSFR) ips = ip/M^, for the AKARI- 
SDSS sample (see also Takeuchi et al. 2010). The SDSS 
galaxies without AKARI detection are bimodally dis- 
tributed in SSFR refiecting the well known populations 
of the red sequence and the blue cloud (e.g., Strateva et 
al. 2001; Bell et al. 2004). The AKARI-detected galaxies 
are distributed in the high SSFR and high ry region, in 
agreement with the general expectation that recent star 
formation and a large dust optical depth are necessary for 
a galaxy to be bright in far-infrared wavelength. However, 
our sample covers a wide range of SSFR and tv, indicat- 
ing that it includes galaxies having a wide range of star 
formation history and dust opacity. 

Finally Figure 10 shows A/ri versus redshift for 
the AKARI-SDSS and AKARI-HIPASS samples. The 
HIPASS catalog is limited by H I flux, and this can clearly 
be seen in the z-Mhi relation in the right panel. The ef- 
fect of H I fiux sensitivity limit is also seen for the AKARI- 
SDSS sample, though it is less clear because of the het- 
erogeneous flux limits of the Hi data in the literature. 
Galaxies detected by AKARI have higher Hi mass on 
average at a fixed redshift than those without AKARI- 
detection, and hence the AKARI sensitivity limit is also 
important for these samples. 
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Fig. 4. Rcdshift versus AKARI band flux for the AKARI-SDSS (left) and AKARI-HIPASS (right) samples. The dashed 

lines indicate the flux sensitivity limit (the flux where the detection completeness is 50%). 




Fig. 5. Angular size versus optical magnitude for the AKARI-SDSS (left) and AKARI-HIPASS (right) samples. The grey dots 
arc galaxies in the original SDSS or HIPASS catalog, without requiring detection by AKARI. The grey dots in the left panel are 
randomly selected from the whole SDSS catalog because of the large sample size of SDSS. 



3. Fitting Physical Dust Models to Infrared 
SEDs 

3.1. The Physical Dust Models 

We mainly use the physical dust model of Draine & Li 
(2007, hereafter DL07) to derive physical quantities about 
dust emission (see also Weingartner & Draine 2001; Li & 
Draine 2001). To reduce the number of model parameters, 
we choose the jm = 5 model of DL07 for the Milky Way 
(MW) that has the PAH mass fraction of gpAH = 3.19%, 
because this PAH fraction is close to the median of nearby 
galaxies studied by Draine et al. (2007). The normaliza- 
tion of the DL07 model is determined by the hydrogen 
mass Afjjj assuming the dust-to-hydrogen mass ratio of 
fdH = Md/Mii ~ 0.01, but fdH likely depends on proper- 
ties of galaxies, especially metallicity. Therefore in this 
work we quote the dust mass as the normalization param- 
eter in the SED fitting, rather than the hydrogen mass. 



The dust mass inferred from infrared SED fittings is de- 
noted by MdiK, to make a distinction from the total dust 
mass Md in a galaxy potentially including cold dust that 
is not heated by radiation and hence not contributing to 
infrared emission. 

The basic parameter of infrared emission in the physi- 
cal dust model is the heating radiation field strength, t/, 
which is the dimensionless energy density of dust-heating 
radiation normalized by the local value of the interstellar 
radiation field (ISRF) around the solar neighbourhood. In 
reality, the field strength can be different at different loca- 
tions in a galaxy, and model infrared emission is often cal- 
culated considering the dust mass distribution dMdm/dU 
as a function of U (Dale et al. 2001; Dale & Helou 2002; 
DL07). However, the functional form of dAIdm/dU is 
highly uncertain, making it difficult to interpret the re- 
sults of infrared SED fittings in a model-independent way. 

Here we take a different approach by assuming that in- 
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Fig. 7. The same as Fig. 5, but for optical color versus absolute optical magnitude. Colors were calculated from absolute magnitudes, 
which have been corrected for the Galactic extinction. The X-eorrection has been done for the SDSS magnitudes but not for the 
HIPASS magnitudes. 




redshift z redshift z 



Fig. 10. Neutral hydrogen mass versus redshift for the AKARI-SDSS (left) and AKARI-HIPASS (right) samples. The grey dots 
are galaxies in the original SDSS or HIPASS catalog, without requiring detection by AKARI. 



frared emission from a galaxy can be described by a single 
characteristic heating radiation strength, Vh- Although 
it may be difficult to fit observed infrared SEDs in a 
wide wavelength range from near- to far-infrared by a sin- 
gle value of [/, in this work we concentrate on the fit- 
ting in a relatively narrow range of wavelength covered 
by AKARI/FIS (65-160 /im), corresponding to the broad 
SED peak of modified blackbody. Therefore, Uh in this 
work should be regarded as the characteristic value for 
large grains contributing to the thermal emission. An ad- 
vantage of this approach is that physical discussion based 
on the radiation field strength becomes easier. The SEDs 
of the single-?/ DL07 model for several values of Uh are 
plotted in Fig. 11. 

It would be convenient to compare the SEDs of the 
physical dust model with the modified blackbody spec- 
tra often used in the literature. It is known that SEDs 
around the thermal peak can be approximated by a mod- 



ified blackbody^ with the emissivity index of /3 ~ 1-2, i.e., 
flux density per unit frequency Fi,(y) cx Qcmiy) B^lv, Td), 
where B^, is the blackbody spectrum of temperature Td 
and Qom oc is the emissivity of dust particles (e.g. 
Dunne et al. 2000; Kovacs et al. 2006). The total 
infrared luminosity of the modified blackbody scales as 
-^TiR oc T^'''^ , and the vF^ peak scales as J^poak ocT^. Since 

-^^TiR <y^Uh for a fixed dust mass, we expect C//i oc t'pjlak- 
The relation between Ut and v-peak of the DL07 model 
implies (3 = 1.75, and for this value the modified black- 
body with Td = 15, 22, 33, and 49 K has the same :^pcak 
as the single-?/ DL07 model of Uh = 0.1, 1, 10, and 100, 
respectively. These spectra arc also plotted in Fig. 11. 

* In the astronomical literature, this is also often referred to as 
"greybody" , but we avoid this word because strictly it is defined 
as emission with a constant emissivity less than one against 
wavelength. 
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Fig. 9. The same as Fig. 8, but for V-band optical depth 
of dust in galaxies (tv) versus specific star formation rate 



Fig. 11. The infrared SED models of galaxies. The red thick 
solid curves arc the single-C/ predictions by the physical dust 
model of DL07, for several values of Ui^ indicated in the fig- 
ure (from bottom to top with increasing (7^). The blue thin 
curves are the same but predictions including the PDR com- 
ponent with 7 = 0.04, for several values of Uaiin- The dotted 
curves are the corresponding PDR components. The dashed 
curves are predictions by the DH02 model for the four val- 
ues of odh (from top to bottom with increasing qdh)- The 
dot-dashed curves are the /3 = 1.75 modified blackbody spec- 
tra for several values of dust temperature (from bottom to 
top with increasing T^). 

3.2. Fitting to the AKARI/FIS Data 

Because the number of galaxies detected by IRC is small 
and we concentrate on Ut for the broad thermal peak in 
SEDs, we use only the four AKARI/FIS bands in the SED 
fits. For a galaxy whose energy flux per unit frequency 
is F^{v), the i-th FIS band flux quoted in the FIS-BSC 
catalog (in units of Jy) is calculated as 



F, 



1 



dvF,{y)R,{y) , 



(5) 



where v is the frequency in the observer's frame, Ai/^ is 
the band width 



Am, I dvv ^ Riiy) , 



(6) 



Vc~c/\c is the central frequency of the band given in 
Table 1, and Ri{v) is RSRF against energy flux for each 
AKARI band shown in Fig. 1. The band width is defined 
so that Fi~F^{vc) when vFy{v) is constant. The observed 
energy flux is related to luminosity Ly{vr] AfdiR, Uh) of 
the DL07 model per unit rest-frame frequency at a redshift 
z as 



(1 + ^) 



Ly{vr; Mdin, Uh) 



(7) 
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where the rest-frame frequency is z/^ = (1 + z) v. 

Then the best-fit values of Mam and Uh are calculated 
by the standard minimization, i.e., 



X 



f: 



'obs 



F- 



model 



(MdIR,f//0]' 



tot,i ) 



(8) 



where the summation is over the four FIS bands. When 
the flux measurement is not reported for a FIS band in 
the FIS-BSC catalog, this band was not used in the fit- 
ting. The number of the FIS bands used for a fitting to 
a galaxy is denoted by A'dat , and because of our require- 
ment of FQUAL > 3 at least in two of the four FIS bands, 
2 < -/Vdat < 4. There are 3, 44, and 831 galaxies having 
N^e^t = 2-4 for the AKARI-SDSS sample and 4, 33, and 
674 for the AKARI-HIPASS sample, respectively. The 
range oiUh covered in the fittings is 0.1 < J7;i < 10^, where 
the lower bound is limited by the available model library, 
and the upper bound is large enough because we found no 
galaxies having Uh > 100. There are 20 and 24 galaxies 
whose best-fit Uh is 0.1 (the lower bound of the model) 
for the AKARI-SDSS and AKARI-HIPASS samples, re- 
spectively, and the following statistical analyses are not 
significantly affected by these small number of galaxies. 
The total infrared luminosity Xtir is then calculated as 
the bolometric luminosity of the dust emission from the 
best-fit SED model. 

The 1(7 error AUh for Uh was calculated for each galaxy 
by the range corresponding to Ax^ = 1 from the min- 
imum, taking the best-fit normalization (i.e., MdiB.) for 
each value of Uh- The errors of M^m and ixiR are cor- 
related with Uh, and we estimated them as follows. For 
a fixed value of Uh-, the error for the normalization factor 
can be estimated as (see, e.g.. Press et al. 2007) 



AM 



dIR 
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1 /I a^x^ 



2 dMj,^ 



In addition, to take into account the error correlated with 
Uh, we calculated the change of best-fit Mdm and Ltir 
values when Uh is changed from the best-fit value by 
±AUh- Then the quadratic sum of these two types of er- 
rors was calculated as the final error for M^ir and Ltir,. 
We found that the mean errors of log^Q Uh and log^Q MdiR, 
are - 0.17 both for the AKARI-SDSS and AKARI- 
HIPASS samples, and that for logj^Q Ltir is smaller than 
these by a factor of about four. The error of Ltir is small 
because it is directly connected to the observed fluxes, 
while MdiR and Uh have similar fractional uncertainties 
about the SED shape by the relation Ltir oc A'/^ir Uh- 

Six examples of the fittings to the AKARI-SDSS galax- 
ies are shown in Fig. 12. The distributions of Uh and 
the minimum reduced chi-square, Xred = X^/'^ ^oi Adat > 3 
galaxies are shown in Fig. 13, where the degree of freedom 
is 1/ = N^nt — 2 for the two-parameter fittings. Theoretical 
expectations of the Xred distribution (weighted sum of 
those for iVdat = 3 and 4) arc also shown. The observed 
Xred distribution is larger than the ideal theoretical pre- 
diction by a factor of 1.5, but we consider that this is 
not unreasonable, if we take into account the possible sys- 
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Fig. 13. (Top) Dust-heating radiation field strength (Uh) 
versus the minimum reduced chi-square (x^cd) °^ SED fit- 
tings for the AKARI-SDSS and AKARI-HIPASS samples. 
(Bottom) The histogram of the minimum X^cd' '^^^ curves 
are the theoretically expected distribution from the statis- 
tics for ly = (A^dat ~ 2) degrees of freedom. Only galaxies with 
-'^dat ^ 3 are shown in this figure. 



tematic errors in the AKARI fiux calibration and/or the 
possible difference of the single-?/ DL07 model from the 
reality. 



(9) 3-3. Models versus Observations in the IRC Bands 



In this paper wc will concentrate on the SED param- 
eter Uh determined by fits to the AKARI/FIS bands, 
but it would also be interesting to compare the observed 
AKARI/IRC fluxes with those predicted by the physical 
dust models. This is done in Fig. 14, where we show the 
flux ratios of FigZ-Fgo and Fg/Fig against F^q/Fi^q, where 
Fj. denotes the AKARI photometric flux of the x-^m 
band. Here we have used galaxies with the AKARI/FIS 
flux quality flag FQUAL > 3 both in the 90- and 140- 
^m bands from the AKARI-SDSS and AKARI-HIPASS 
samples. The ratio of i^9o/^'i4o can be regarded as an in- 
dicator of the dust temperature or Uh- The solid curves 
are the path along Uh predicted by the single-C/ DL07 
model, and the prediction for Fi^/Fgo is lower than the 
mean of the observed distribution by a factor of about two. 
This is reasonable because mid-infrared emission is likely 
dominated by dust grains in regions of stronger heating 
radiation fleld, and we must incorporate the diversity of 
U to fit SEDs in a wide wavelength range. Here we adopt 
two such SED models. 

One is the physical dust model by Dale & Helou (2002, 
hereafter DH02), which introduces a power-law distribu- 
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Fig. 12. Examples of the SED fittings to the galaxies in the AKARI-SDSS sample. The squares are AKARI data points, and error 
bars arc small and difficult to see for most of them. The three different models are used for the fits: the 7 = DL07 model with the 
parameter Uh (solid), the 7 = 0.04 DL07 model with the parameter C/min (dashed), and the DH02 model with the parameter odh 
(dotted). The best- fit parameters and X^ej = X^/'^ (with 1/ = N^^^ — 2 degrees of freedom, for the 7 = DL07 model) are shown in 
each panel. The unique identification set of plate number, MJD, and fiber number of the SDSS spectroscopic database is also shown 
in each panel. 



tion of radiation strength as 



dU 



(10) 



and the path along odh predicted by this model is shown 
in Fig. 14. The values of aDH = 4,3,1.8, and 1 roughly 
correspond to Uh = 0.6,1.0,10, and 100 of the single-[/ 
DL07 model, respectively, in our SED fits to the four 
FIS bands. This model shows a better agreement with 
the data than the single-t/ DL07 model, demonstrating 
the necessity of including a small amount of large U dust 
component. On the other hand, the data show dispersions 
by about one order of magnitude in both the flux ratios 
of Fg/Fis and Fis/Fgo, indicating the limitation of one- 
parameter-family SED models. The SEDs of the DH02 
model arc also shown in Figs. 11 and 12. 

The other model is DL07 using a form of dMdin/dU 



proposed by DL07, which introduces a power-law distri- 
bution of dAIdm/dU in the range of J/min <U < J7max for 
a fraction 7 of the total dust mass, in addition to the 
single-?/ component with U — Umin, i.e., 



dM, 



dU 



= (l-7)MdiR5((7-[/,„i„) 



a-l 



dlK 



^ max 



u- 



(11) 



The power-law component is called as the PDR compo- 
nent as it represents the hot dust in photodissociation 
regions (PDRs) close to young massive stars. We adopt 
a = 2 and J/max = 10^ for the PDR component following 
Draine et al. (2007), and now we have another model pa- 
rameter 7 in addition to C/min- (The model has yet another 
degree of freedom about the PAH fraction by the model 
parameter jm , but we keep Jm = 5 to reduce the param- 
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eter space to explore.) The PDR component models are 
available only for t/min < 25, and we use the Umm = 25 
template for the PDR component when C/min > 25. We 
calculated the path along Umin in the plane of Fqq/Fuo 
versus Fg/Fig and Fig/FgQ for various values of 7, and 
found that 7 = 0.04 gives the median of the observed dis- 
tribution, as shown in Fig. 14. This value is close to 
the typical value found by Drain et al. (2007) for nearby 
galaxies, though there is a significant scatter of 7. The 
SED model curves of the 7 = 0.04 model are also shown 
in Figs. 11 and 12. 

Although the 7 = 0.04 model is better than the single-?/ 
model (i.e., 7 = 0) to reproduce the observed IRC fluxes, 
we will use the latter in the analyses below because of 
the following reasons. First of all, we want to estimate 
the characteristic radiation field strength Uh responsible 
for the emission around the thermal peak from large dust 
grains, and introducing the PDR component in a model- 
dependent way may induce some bias in the estimates of 
Uh- The 7 = 0.04 model produces a feature around 30 fim 
because of the rather arbitrary parameter of C/max = 10^, 
when Uaiin is small. The best-fit Umin values of the 7 = 
0.04 model tend to be smaller than Uh of the 7 = model 
(see Fig. 12), and we found that a considerable number 
of galaxies have the lowest value of Umin = 0.1 provided 
in the model library. The number of such galaxies can be 
minimized by using the 7 = model. Finally, it should 
also be noted that the deficit of the predicted IRC-band 
flux by the 7 = model is at most a factor of two compared 
with the mean observed fiux. 

4. Theoretical Considerations 

Before we compare the infrared SED parameters with 
other physical quantities of galaxies, we describe the gen- 
eral theoretical expectations. The basic assumptions here 
arc that dust and heating radiation sources are uniformly 
distributed on a disk (but possibly with different scale 
heights). This picture is obviously simple, and allowed 
only when the infrared SED is determined mostly by the 
general ISRF on the global scale in a galaxy, rather than 
those associated with small scale structures like individ- 
ual star forming regions. However, rather surprisingly, we 
will find that this simple picture gives reasonable expla- 
nations for most of the data presented in the paper, and 
more geometrically complex models are not warranted. 

Basic Equations of Emission from Dust 

The infrared SED depends on the heating radiation field 
strength felt by dust particles, which is defined by the en- 
ergy density u,^ or energy fiux Si/ per unit frequency, and 
these two are simply related where c is the 

speed of light. The radiation field can be anisotropic, 
but anisotropy does not affect the infrared emission from 
heated dust, and these quantities are defined as those inte- 
grated over 47r steradians. In this work we do not consider 
the variation of the SED of the heating radiation field, and 
use the same parameter Uh introduced in §3, which is a di- 
mensionless field strength normalized by the local strength 
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Fig. 14. The color index correlations for the AKARI-SDSS 
and AKARI-HIPASS samples. The flux denotes photo- 
metric flux in the AKARI x-^ra band in units of Jy. The 
three curves are the predicted path along the model parame- 
ters by the physical dust models. The solid, dashed, and dot- 
dashed curves are for the 7 = DL07 model (Uh = 0.1-100), 
the 7 = 0.04 DL07 model = 0.1-100), and the DH02 

model (ojjH = 0.0625-4), respectively. The model predictions 
assume z = 0.02 that is typical for the observed galaxies. 
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u^fi of the Milky Way ISRF around the solar neighbour- 
hood, i.e., u^{i') = Uh u^fliu) and Si^iu) = Uh Syfiiu). 

When dust particles of total mass Af^m are exposed to 
a heating radiation field strength Uu and emitting in a sta- 
tionary state (a good approximation in the AKARI/FIS 
bands), the bolometric luminosity Ltir emitted from 
these particles is given by the energy balance with the 
total energy absorbed by dust grains, as 



^TTR — M, 



dm. 



dv Kd{v) Sy{v) 



(12) 
(13) 



where is the dust mass opacity coefficient (absorp- 
tion cross section per unit mass of dust grains, see, e.g., 
Hildebrand 1983; Dunne et al. 2003). If we define the 
frequency-integrated effective dust mass opacity weighted 
by Si,fi, as 



J dvs„fi{v) 



(14) 



we obtain 



Ltik = Kd,c« MdmUh s\,o\,o , (15) 

where Sboi.o = / s^.o dv is the bolometric intensity of the 
local ISRF, for which we adopt the value of Mathis et al. 
(1983): Sboho = 0.0259 erg cm-^s"! in 0-13.6 eV (see also 
Draine et al. 2007). We can calculate K^^cff for the physi- 
cal dust model of DL07 from eq. (15), which is found to 
be Kd^cS = 1.0 X 10^ cm^g^^. In astronomically convenient 
units, this can also be written as 



^^1.4xl02^C/, 



(16) 



Now, adopting the disk geometry, we define the global 
infrared radiation field strength of a galaxy by 



1 -^TIR 1 

Stir, = x — ~ 7: nd.cs ^dm Uh Sboi,o 



(17) 



where R is the characteristic disk radius of a galaxy, and 
the dust mass column density of a galactic disk is de- 
fined as EdiR = Mdm/ {-kR^). In this work we use the 90% 
Petrosian radius Rpgo as R for the AKARI-SDSS and 
AKARI-HIPASS samples, which includes most of the to- 
tal luminosity from a galaxy. For the AKARI-SB galaxies, 
we use the size of circumnuclear starburst regions given 
in K98. The factor of 1/2 is introduced to approximately 
account for the disk geometry where the radiation is emit- 
ted into the two sides of the disk. The quantity Etir has 
the same physical dimension as the heating radiation field 
Sboi = Uh Sboi,Oj and we expect some relation between the 
two if the infrared emission is mainly determined by the 
global properties of a galaxy. To directly compare with 
the dimensionless strength Uh, we normalize Etir as 



f^TIR : 



^TIR 
Sbol.O 



1 



Kd,cff ErfIR Uh 



(18) 



The effective optical depth Td.cS of dust, corresponding 
to the light travel distance I, is defined as r^^eff = Kd,eSPdl, 



where pd is the spatial mass density of dust in ISM. The 
optical depth is related to dust mass column density and 
V band extinction as: 




= 0.327 



(19) 



where = 0.48 Mq pc ^ and we have used the relation 
between Ay and hydrogen column density, Av/Nyi = 5.3 x 
10~^^ mag cm^ (H atom)~^, for the DL07 physical dust 
model. We also introduce the total optical depth t^I^ of 
the disk, which is defined as 



_tot ^ ^ 



1 Kd^cfiMd ^ _ 

— 7; ^d,cff ^d — 



1 

2 ' 



E(i,crit 



(20) 



where Md is the total dust mass in a galaxy, 'Sd = 
Md/ [ttR^), and we have defined the critical dust mass 
column density E^.crit = 2 n^^a ^ ^-^^ Mq pc^^ for which 
T^°^ff = 1. The factor 1/2 has been introduced because, 
on average, radiation from heating sources passes through 
about half of the total column density of the dust disk. It 
should be noted that here we make a distinction between 
Mdm (heated dust mass estimated from infrared emission) 
and the total dust amount Md in a galaxy. Generally 
-^diR (or E^ir) can be smaller than Md (or E^), if only a 
fraction of dust particles are radiated by strong radiation 
field. The total disk optical depth T^°*ff is defined by Md 
rather than Mdm- Then we have 



C^TIR - 



E(JIR 



rr _ tot Mdm. 



Uh 



(21) 



Ed '•^^ Md 

This is an important equation that will be used to inter- 
pret the observed data in this work. 

4-2. Settlings among [/tir, Uh, Mdm, o,nd Aid 

Eq. (21) is a generic equation, but we can consider some 
limiting cases with simple scaling laws, corresponding to 
a particular geometrical distribution of dust and heating 
sources. The schematic picture to explain these cases is 
given in Fig. 15. 

First we consider the case of small optical depth, T^°^ff ^ 
1 (hereafter the case A). In this case the heating radia- 
tion field strength should be Sboi — Uh Sboi.o ^ Uh/ {2nR'^), 
where Lh is the bolometric luminosity of heating radiation 
from all sources in a galaxy. This field strength is relevant 
for all dust particles in a galaxy, and hence M^ir ~ Md, 
while only a small fraction of the heating radiation energy 
is absorbed and converted into the infrared emission, i.e., 

t^TIR ~ tJ^cS < Uh- 

In the opposite limit of large optical depth (T^°*ff ^1), 
most of heating radiation will be converted into infrared 
emission, and hence Ltir ^ Lh and Etir ~ Lh/ {2ttR^) 
(hereafter the case B). We can further consider two very 
different limiting cases about the relative distributions of 
dust and heating sources. The first case is when dust 
particles have a similar disk scale height as the heating 
sources. Then all the dust particles in the disk contribute 
in the same way to the infrared emission, i.e., M^ir = Md- 
Then eq. (21) tells us {/tir = rf *ff Uh, i.e., JJtir > Uh, 
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Fig. 15. Schematic pictures for the representative cases of the geometrical distributions of dust-hcating stars and dust grains, in 
the disk geometry. The normal of the disk plane is along the vertical direction. 



meaning that the heating radiation strength felt by dust 
particles (Uh) is smaller than Utir = Lh/ {^-kE? ShoXfi) by 
a factor of T^°eff. This is because the heating photons are 
substantially absorbed after traveling through a distance 
corresponding to Td,cS — l, and only a fraction (~ l/'''rf°off) 
of total heating radiation produced in the disk reaches a 
dust particle. We call this as the case Bl. 

Another optically thick case is when the disk scale 
height of heating sources is smaller than that of dust parti- 
cles. As a simple example, we assume that the scale height 
of heating sources is negligibly small, so that the dust disk 
behaves like an optically thick screen. Then the radiation 
field strength at regions close to the heating sources should 
be Sboi = Lh/i^TiR^) = Etir, i.e., Uh = Utyr. Then eq. 
(21) requires 

meaning that only a fraction '^/t^cs) '^^ViSt particles 
are heated by the radiation, because the radiation field 
is significantly damped after traveling through a distance 
corresponding to Td.cS = 1- Infrared emission is mostly 
from the dust in the thin layer corresponding to r^j cff ^ 1 
around the heating sources. Then it is predicted that the 
column density of heated dust is constant at I]c(ir~ Sd^crit, 
regardless of the total dust column density S^. We call 
this as the case B2. 

Note that M^ir ~ Md and [/tir ^ Uh can be realized 
simultaneously only when T^°*ff ~ 1- 

5. Results and Interpretation 

Here we present the main results of this work, i.e., the 
correlation between infrared SEDs and various physical 



properties of galaxies. The results will be interpreted by 
the theoretical background described in the previous sec- 
tion. 

5.1. Infrared Luminosity versus Optical Luminosity and 
SFR 

First we compare the infrared luminosity Ltir with 
the optical luminosity, to examine the dust opacity T^°*ff- 
Figure 16 shows the correlation between Ltir and the op- 
tical r-band luminosity {Lr = vL^, at rest-frame 6231 A). 
We find Ltir ^ Lr for most galaxies, indicating that en- 
ergy output by dust emission is similar to, or larger than 
that by the direct stellar emission. It may indicate that 
most galaxies have r^°efj ^ 1, but we can derive this conclu- 
sion only if the LxiR/^r ratio is determined by the global 
properties of galaxies, because r^°*g has been defined as 
the mean dust column density on a galactic scale. Instead, 
the Lj^m/Lr ratio may be determined by the properties 
of local star forming regions in a galaxy. We will find that 
the former interpretation, rjf^ff ^ 1, is indeed appropriate 
for the galaxies studied here. 

Next we examine the correlation between Ltir and SFR 
(V', estimated from SDSS spectra) of the AKARTSDSS 
sample, shown in the left panel of Figure 17, to examine 
the origin of dust heating radiation. The total infrared 
luminosity Ltir is often used as a SFR indicator by eq. 
(4), and its validity and limitation have been discussed 
in the literature (e.g., Takeuchi et al. 2005, Salim et al. 
2009; Boquien et al. 2010; Goto et al. 2011b). It cannot 
be a good SFR indicator when star formation activity is 
low and dust heating radiation is significantly contributed 
from relatively aged stars. This trend is clearly seen in 
this figure; Ltir/^ increases with decreasing SFR. The 
solid line is the Xtir/^ ratio of eq. (4), which is the 
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value expected when the dust heating radiation is dom- 
inated by young stellar populations with ages of 10-100 
Myr (Kennicutt 1998, converted to the Kroupa IMF here). 
This Ltir/'0 ratio is observed only for high SFR galax- 
ies. Different symbols are used in this figure for differ- 
ent ranges of LxiR/ir-j and the trend of larger LxiRZ-^r 
for larger SFR can be seen, in agreement with the well- 
known fact that actively star-forming galaxies are gener- 
ally dusty. 

Since SSFR {ips = tp/M^,) is a good measure of relative 
abundances of young and old stellar populations, we ex- 
pect that Ltir/V' is tightly correlated with SSFR. This 
is indeed seen in the right panel of Figure 17, especially 
for galaxies having L^m/Lr < 10. On the other hand, 
dusty galaxies having L'yyrJ Lr > 10 deviate from the cor- 
relation showing particularly high Ltir, and this is most 
likely explained by the hidden activity of star formation 
that cannot be captured by optical emission because of 
large extinction. The tight correlation of galaxies with 
Ltir/Lt < 10 can be fit as 

This result indicates that SSFR can be used as an in- 
dicator of contribution from young stellar population to 
dust-heating radiation. The infrared luminosity Xtir can 
be used as a simple SFR indicator only for actively star- 
forming galaxies of -08 ^ 0.1 Gyr~^, while for other galaxies 
relatively aged stars significantly contribute to the dust- 
heating radiation. The tight Xtir/V' versus ips correla- 
tion may be useful for parameter estimations for observed 
galaxies, because if two of the three parameters (M*, SFR, 
and ixiR.) are known, the remaining parameter can be de- 
duced by solving the Ltir/'0 - "^s relation. For example, 
if one knows Xtir and Af* from observations, SFR can be 
solved from this relation, which should be a better esti- 
mation than simply using eq. (4). 

5.2. Infrared SED versus Infrared Luminosity and SSFR 

Figure 18 shows the correlation between Ltir and Uh, 
for the AKARI-SDSS and AKARI-HIPASS samples. As 
mentioned in §1, the correlation between the two (larger 
Uh for larger Ltir) is widely known, and the trend is 
also confirmed for both the two samples studied here. 
However, there is a considerable scatter and it is difficult 
to predict Uh with a good accuracy only by the informa- 
tion of Ltir. Furthermore, we find a significant offset 
of the mean relation between the two sample. Since the 
AKARI-SDSS galaxies are more distant, their mean Ltir 
is larger than that of the AKARI-HIPASS galaxies, but 
the distribution of Uh is not significantly different. This 
indicates that the ixiR-f^/i or LTiR-7d relation is sensitive 
to the sample selections. 

Figure 19 shows the correlation between Ltir and aoH 
of the DH02 model fits. It is interesting to compare the 
aoH distribution of our samples with that of the IRAS 
Fqq > 1.2 Jy sample (Fisher et al. 1995). Converting the 
distribution of Fqq/Fiqq flux ratio of this sample given 
in Chapman et al. (2003) into q;dh using the relation 




specific SFR [Gyr 



Fig. 20. Dust-heating radiation strength (Uh, estimated 
by infrared SED fittings) versus specific star formation 
rate (ips = t/j/A'It, iji and M» from SDSS spectra) for the 
AKARI-SDSS sample. 

in Table 2 of Dale et al. (2001)^, the mean is odh ^ 2 
and there are few galaxies having aoH > 2.5 in this IRAS- 
based sample. This is in contrast to the AKARI-SDSS 
and AKARI-HIPASS samples that include many galaxies 
at aDH > 2.5. This is likely because the IRAS sample 
was selected at the IRAS 60-/im band, while many of the 
galaxies in our sample are detected in the AKARI 90- 
and 140-/Ltm bands with the FIS flux quality flag EQUAL 
> 3 (see Table 3). Therefore our samples include galax- 
ies having colder dust temperatures than the IRAS-based 
sample. Note that the errors of q;dh (calculated in the 
same way as Uh) are large when q;dh ^ 3, because the de- 
pendence on ckDH of the DH02 templates becomes weaker. 

Figure 20 shows the correlation between SSFR and Uh 
for the AKARI-SDSS sample, using SDSS SFR and stel- 
lar mass. This correlation is physically more reasonable 
than LxiR-t^h, because both SSFR and Uh arc intensive 
physical quantities (i.e., not depending on the system size 
or the amount of material). In fact, we see a better corre- 
lation in Fig. 20 than the LTiR-Uh correlation: larger Uh 
for galaxies having larger SSFR. However the slope of the 
correlation becomes steeper at tps ^0.1 Gyr~^, and it is 
still difficult to predict Uh reliably simply from SSFR. It 
should also be noted that SSFR is not a directly observ- 
able quantity, and physical relation between SSFR and Uh 
is not clear. SSFR is related to star formation history in a 
galaxy and hence it would affect the SED of dust-heating 
radiation field. However, a key quantity controlling the 
radiation field strength is the geometrical size, which is 
not directly related with SSFR. Correlations between Uh 
and more directly observable quantities based on better 
physical connections are worth to seek for. 

^ Here we ignored the X-correction, which is a good approxima- 
tion for the low-redshift sample. 
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star formation rate if/ [M^ yr '] specific SFR ij/^ [Gyr '] 



Fig. 17. Correlation between SFR (tp, estimated by SDSS spectra) and Ltir (left), and that between SSFR (ips ='4>/Mt) and 
^tir/V' (right) for the AKARI-SDSS sample. Different symbols are used for different ranges of the infrared-to-optical luminosity 
ratio, Lq^i^/Lr, as shown in the left panel. The solid lines are the calibration of Ltir/V" by Kennicutt (1998, but corrected for 
the Kroupa IMF) for -Ltir as a^n SFR indicator. The dashed line in the right panel is the best-fit power-law to the galaxies with 
LTm/Lr < 10. 
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Fig. 18. Dust-hcating radiation field strength (Uh, estimated by infrared SED fittings) versus total infrared luminosity LtiR; for 
the AKARI-SDSS (left) and AKARI-HIPASS (right) samples. 




total IR luminosity log^g L^^^^ [Lg] total IR luminosity logjg L^jg [Lg] 

Fig. 19. The same as Fig. 18, but for the correlation of qdh (obtained by the SED fittings with the DH02 model) versus Ltir- 
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5.3. Infrared SED versus Global Infrared Radiation Field 
Strength 

Figure 21 shows the correlation between Uh and ?7tir, 
i.e., the galactic-scale infrared radiation field strength 
Etir = LTiR/(27ri?^) normalized by the solar neighbour- 
hood value of ISRF. It can be seen that these two pa- 
rameters are correlated much better than the L^iRrUh or 
ips-Uh relations. The relation is well described by the lin- 
ear relation of Uh ex Utir, with the ratio Uh/UTiR, ^ 1. 
This result is in nice agreement with the prediction in the 
case B2 discussed in §4.2, and it strongly indicates that 
(1) integrated infrared SED of a galaxy is mostly deter- 
mined by the general ISRF on the global scale of a galaxy, 
and (2) dust-heating sources are deeply embedded in an 
optically-thick dust disk, i.e., the scale height of heating 
sources being smaller than one optical depth of the dust 
disk. 

There are several independent observational evidences 
supporting this picture. Walterbos & Greenawalt (1996) 
and Sun & Hirashita (2010) argued for the dominant role 
of the general ISRF in the determination of infrared SEDs. 
The geometry of the case B2 is consistent with the smaller 
scale height of molecular gas than that of H I in the Galaxy 
(Malhotra 1994, 1995; Binney & Merrifield 1998), consid- 
ering that young stars are generally born in dense molec- 
ular gas. Here we implicitly assumed that dust-to-gas 
ratio is fairly constant throughout the Galaxy, which is 
reasonable as inferred from the well-known good corre- 
lation between total (i.e.. Hi plus H2) hydrogen column 
density and color excess of extinction (Bohlin et al. 1978; 
Heithausen&Meboldl989;Rachfordetal. 2009). Finally 
we note that the original form of the Schmidt law (?/; cx p^^^ 
where pgas is spatial gas density) was proposed to explain 
the smaller vertical scale height of young stars than that 
of interstellar gas (Schmidt 1959; Fuchs et al. 2009). 

In the case B2, we expect that the dust column density 
inferred from infrared SED fits, S^ir, is constant around 
Sd.crit = 2 K-^l^ = 0.96 [Mq pc-2]. This is checked by Fig. 
22, where correlation between C/tir and E^ir is shown. 
As expected, the values of E^ir do not show any trend 
against C/tir, with the mean value close to E^^crit- Many 
of galaxies having large deviations from Ed,crit have large 
observational errors. Note that no galaxies are found in 
the region of ?7/i < 0.1 simply because the fit is limited by 
the minimum value Uh = 0.1 of the physical dust model. 

We fit the C/tir- C//i relation by a linear {Uh = 10° C/tir) 
relation and the results are presented in Table 4 for the 
AKARI-SDSS and AKARI-HIPASS samples. Here we 
considered only the errors of logi^Uh, because the errors 
of logj^o^TiR are much smaller, as mentioned in §3.2. The 
proportionality constant is different by just a factor of 1.5 
for the two samples, indicating that our result is not seri- 
ously affected by the selection bias. We also tried a power- 
law {Uh = IO^CZ-^ir) relation and a three-dimensional plane 
fit in the space of logj^pi?, logj^piTiR, and logjg C//i. The 
results are shown in Table 4, and these fits also indicate 
that the simple relation of Uh oc C/tir oc L^ir/B? is a good 
description of the observed correlation. 



Table 4. Functional Fits to the C/TlR-f/fe Relation 
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We estimated the physical dispersion of Uh from the 
mean relation as follows. We define A°^^ as the observed 
deviation of log^QUh of each galaxy from the best-fit rela- 
tions. Then the observed median of |A°]^^| was converted 
into the observed standard deviation cr°g^^ by using the 
Gaussian relation (median of \x\ is 0.674 cr). The use of 
the median is to avoid the effect of a small number of 
outliers having large |A°^^|. The quantity (J°^lr includes 
the observational errors of Uh in the SED fittings, and 
the intrinsic (i.e., physical) standard deviation crjg[/ was 
calculated as 

<c/=[«a)'-(4V)']' , (24) 
where (cf^V) ^^'^ mean of la errors of log^o Uh in the 
SED fits. The resuhs for cr°^l,, a{f,j and (af^V) are shown 
in Table 4, and o-"gi; ^0.3 is indicated for all cases. This 
corresponds to ^13 % dispersion in terms of the modified 
blackbody temperature, using Uh^xT^^^ and /3 = 1.75. 

5.4- Dust Mass versus Metal Mass 

If the interpretation of the case B2 distribution is cor- 
rect for the observed Urm-Uh correlation, this can fur- 
ther be tested by checking the trend about M^ir against 
Md, since the case B2 predicts the scaling of Mdm./Md ~ 
(r^°*jj)~^ oc E^^. For this purpose we need to measure the 
total dust amount in a galaxy independently of the in- 
frared emission. Though it is difficult to directly measure 
Md, the total metal mass of a galaxy in the gas phase, 
Mz, is an indicator of Md because grains are made from 
metals. For a part of the AKARI-SDSS sample, the gas- 
phase metallicity Z and Mm measurements are available, 
and we can estimate the H I-gas-phase metal mass as 

Mzijii) = 1.36 ZMhi = 1.36Z/hiA/h = /hiM^ , (25) 

where /hi is the Hi fraction of the total hydrogen gas 
mass, Mh. [The numerical factor corrects the hydrogen 
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Fig. 21. Dust-heating radiation field strength (Uh, estimated by infrared SED fittings) versus global infrared radiation field strength 
[Etir = -LTlR/(27r_R2)], for the AKARI-SDSS (left) and AKARI-HIPASS (right) samples. The parameter C/tir is Stir normalized 
by the bolometric strength of the local (solar neighborhood) ISRF in the Milky Way, which can directly be compared with Uh- The 
solid and dashed lines are the best-fit linear {Uh oc C/tir) a^nd power-law [U^ oc C/^ir) relations, respectively. 




Fig. 22. Column density of heated dust E^jir versus global infrared radiation field strength {Utir or Sxir for the bottom or top 
ordinate), for the AKARI-SDSS (left) and AKARI-HIPASS (right) samples. The critical dust column density, Ed.crit = 0.98 AfQpC^ 
corresponding to the effective dust optical depth T^°^g = 1, is indicated by the dashed line. The solid line shows the locations when 
Uh = 0.1, and there is no galaxies in the upper-left region because the physical dust model is limited to Uh > 0.1. The errors are 
from those for M^jj^. 
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mass fraction of the total baryon, X = 0.74 (Grevesse et 
al. 2010).] Some of the AKARI-SDSS galaxies and all 
of the AKARI-HIPASS galaxies have Mhi measurements 
but no mctallicity measurements. To calculate Mzini) for 
these galaxies, we adopt a mean value of = 0.05 as in- 
ferred from the mass-metallicity plot of the AKARI-SDSS 
galaxies (Fig. 8). This metallicity is higher than the solar 
abundance, but not unreasonable for dusty galaxies se- 
lected by AKARI, taking also into account the systematic 
uncertainties in the metallicity measurements (see foot- 
note 2). 

Then the total dust mass Md ~ fdz Mz is related to 
A^z(Hi) as 



12 + logjo (0/H) 
8.8 9 9.2 



Md — fdz fm A/z(Hi) 



(26) 



using the dust fraction fdz of the total metal mass in 
interstellar gas. Observations of interstellar medium in 
our Galaxy suggest Md/Mn ^ 0.0073, corresponding to 
fdz ^ 0.3 (Draine et al. 2007). Although /hi is unknown 
for our samples, we expect Md/Mz(Y{i) ^ 1 for galaxies 
with /h I of order unity. 

Figure 23 plots MdiR. against A'/z(hi) for the two sam- 
ples. It can be seen that many galaxies are located at the 
region of Mdm < Mz^m) , and a clear boundary can be seen 
at MdiR Mz(iii), as expected. There are a small num- 
ber of AKARI-SDSS galaxies showing high Mdm/Mz{m}, 
but most of them are galaxies without metallicity esti- 
mates and plotted using Zav. Most galaxies with metal- 
licity measurements have Mdm ^ Mziai)- 

There are galaxies showing Mdm/Mz(m) ^ 1, and the 
origin of this small value is interesting. A possibility is 
that the dust production efficiency from metals is low, 
i.e., Md ^ Mz(Hi)- There are observational indications 
that low-metallicity dwarf galaxies have lower dust-to- 
metal ratios than giant galaxies (Galliano et al. 2003, 
2005; Hunt et al. 2005; but see also James et al. 2002). 
Draine et al. (2007) also found that dust-to-metal ratio 
is lower for low-metallicity galaxies having ^8.1, but 
they argued that infrared emission from such galaxies is 
from a small region compared with the whole galaxy, and 
the dust-to-metal ratio becomes consistent with the MW 
value if it is calculated using gas mass and metallicity in 
the regions of infrared emission. In any case, galaxies in 
our sample have Aq ^ 8.1 as we showed in Fig. 8, and 
hence the low metallicity effect is unlikely to be the origin 
of the small Mdm/Mz^ui) values in our sample. Indeed, 
we plot Mdm/Mz(m) versus metallicity for the AKARI- 
SDSS galaxies with available metallicity measurements in 
Fig. 24, and there is no systematic trend. 

Instead we argue that the dispersion in MdiK/Mz(B_i) 
originates from the difference between Mdm and Md- The 
good correlation between C/tir and Uh indicates the case 
B2 distribution where only a part of dust is radiated by 
heating sources. In this case we expect Mdm/Md oc E^^, 
and hence Mdm./Mz(Hi) « ^z}m) f^z f^l is constant 
for all galaxies, where ^^(hi) = A/2(Hi)/(7ri?^). Figure 
25 shows this correlation, and the data of both sam- 
ples clearly show the expected correlation. This pro- 
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Fig. 24. The ratio of M^i^/Mz(^iii) versus metallicity for 
the AKARI-SDSS galaxies with available Z and Mh I mea- 
surements. The errors are from those for AfdiR- 

vides a further support to the interpretation of C/tir- 
Ufi relation by the distribution case B2. Given the con- 
stant nature of E^iR found in Fig. 22, the correlation 
MdmjMz^Hi) = SdiR/Ez(Hi) « E^(Hi) "^^^^ ^^^^^^ rather 
trivial, but this figure demonstrates that E^iR stays con- 
stant even if E^jhi) significantly changes. According to 
this interpretation, we expect Mdm/Md ^ 1 when T^°^ff ^ 
1, i.e., Tid ~ E£i,crit = 0.98 M© pc~^. The data is consistent 
with this expectation for a reasonable value of fdz /hi^ ^ 1 

(Sz(Hl) ^ ^d)- 

Finally let us examine the assumption of fdz f^^l being a 
constant. The assumption of constant fdz is theoretically 
reasonable if the dust production process works similarly 
in all galaxies. It is also known that the dust-to-gas ra- 
tio estimated by extinction in optical bands for nearby 
galaxies and quasar absorption line systems (QALSs) is 
correlated well with metallicity, indicating that Md/Mz 
is roughly constant (Issa ct al. 1990; Boissier et al. 2004; 
Vladilo et al. 2006). It should be noted that the ex- 
tinction in QALSs is sensitive to all dust grains in a 
system on the line of sight, in contrast to Ad dm that 
traces only dust heated by radiation. The variation of 
/hi in different galaxies may also lead to the scaling of 
MdmjMz(]n) E^^^jjj^. However, if this effect is signifi- 
cant, we expect that the sequence will appear in the region 
of MdmjMz{]i\) ^ 1 because /hi < 1, which is opposite to 
the observed trend. 

5.5. Extension to Circumnuclear Starhursts 

We have seen that the infrared emission from the 
AKARTSDSS and AKARI-HIPASS galaxies are described 
well by the case B2 of dust distribution. However, the dy- 
namic range of U^m is limited to 0.1 ^ J7tir ^ 100, and 
here we examine the AKARI-SB sample to extend this 
range. 
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Fig. 23. Heated dust mass M^jir versus Hl-gas-phase metal mass M^fiu) for the AKARI-SDSS (left) and AKARI-HIPASS (right) 
samples. For red filled triangles in the left panel, Af2(Hi) is calculated using metallicity estimates for each galaxy, while that of open 
symbols is calculated using a mean metallicity of Zav = 0.05 because metallicity measurements arc not available. 




Fig. 25. The ratio of A/diR/Af2{Hi) versus column density of Hl-gas-phase metal ^^{Hi) foi' the AKARI-SDSS (left) and 
AKARI-HIPASS (right) samples. Different symbols have the same meanings as those in Fig. 23. The errors are from those 
for A/diR. 
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Figure 26 shows the Ltir-C/^, Utir-Uh, and J7tir- 
SrfiR relations of the AKARI-SB galaxies in comparison 
with the AKARI-SDSS sample. The AKARI-SB galax- 
ies are distributed at similar locations to the AKARI- 
SDSS galaxies in the L'Ym-Uh plane, but they have much 
larger ?7tir- An interesting trend is found for Uh against 
C^tir; it is roughly constant and there seems an up- 
per boimd of C//i. max ~ 50, which is in sharp contrast to 
the AKARI-SDSS galaxies showing a tight correlation of 
Uh (X C/tir- But these two data sets are not inconsis- 
tent with each other. The dispersion of the Uh/U^iB. 
ratio for the AKARI-SDSS galaxies becomes larger at 
Uh ^ 10, and the distribution of the AKARI-SDSS galaxies 
is smoothly connected to that of the AKARI-SB sample. 
Corresponding to the constant nature of Uh against C/tir, 
SdiR is roughly proportional to C/tir, as expected from 
eq. (18). 

A remarkable result is obtained by comparison between 
Mdin and H2 mass estimated from CO observations, Mh2, 
for the AKARI-SB sample. These two independently mea- 
sured quantities are correlated much better, as shown in 
the histogram of MdiR/MH2 in Fig. 26, than the M^ir- 
Mhi correlation of the AKARI-SDSS galaxies. The ma- 
jority of galaxies are clustering around Mdm/Mii2 ^ 0.015 
within a factor of 2, and the ratio of 0.015 is close to 
the expectation of typical dust-to-hydrogcn ratio fdn — 
1.36fdzZ. It is known that molecular hydrogen is the 
dominant form of hydrogen when the gas column den- 
sity is very high like circumnuclear starbursts (Kennicutt 
1998; Bigiel et al. 2008). Therefore, this resuh implies 
that the infrared-emitting dust mass is close to the total 
dust mass, i.e., M^ir ^ Md- 

We consider that this result can be interpreted as fol- 
lows. As argued by Kennicutt (1998), Ltir can be re- 
garded as a good SFR indicator for the K98 SB sam- 
ple. Using the relations of M^ir ^ 0.015 Mh2 and Xtir cx 
MdinUh, Uh is equivalent to the star formation efficiency 
(SFE), ipe, which is defined as SFR per unit hydrogen gas 
mass. Using eq. (4) and (16), we find 

V'e = = 0.24 Uh [Gyr-i] , (27) 

and the upper limit of C/ft.max ^ 50 is translated into the 
upper limit of V'e.max 10 Gyr^^. Then, the systematic 
change of the C/tir-C//i relation can be explained if there 
is a fundamental upper limit on SFE by the physics of 
star formation (see §6.2 for more discussion on the origin 
of this upper limit). 

If this upper limit exists, the distribution case B2 is not 
allowed when C/tir > Uh.max- For the case B2, star for- 
mation must occur within a thin layer of the disk whose 
dust column density is ^ S^^crit (or Sh Sn.crit 
10^ Mq pc~^, adopting a typical dust-to-hydrogen ra- 
tio). Then there is a maximum of the surface SFR 
density [Ssfr = ''P/{t^R'^)] for the case B2 determined 
by V'e,max and En.crit, corresponding to C/tir ~ Uh,ma,x- 
When SsFR is higher than this, the gas column density 
of the star forming layer must proportionally increase be- 
yond Sn.crit because of the limit of i/'e,max, and hence 



EdiR oc SsFR oc C/tir while Uh becomes constant at Uh.max- 
This situation corresponds to the case Bl rather than B2, 
and the transition from B2 occurs around C/tir ~ CZ/j^max- 
This result should have some implications for the well- 
known Kennicutt-Schmidt law on the plane of Er-Ssfr- 
We will discuss about this issue in more detail in §6.2. 

6. Discussion 

6.1. Paucity of Optically Thin Galaxies and Radiative 
Feedback? 

We have found that galaxies studied here are well de- 
scribed by the distribution case Bl or B2, meaning that 
they are optically thick to dust-heating radiation on a 
galactic scale {rjf^g ^ 1). This is consistent with the dis- 
tribution of the infrared-to-optical luminosity ratio (Fig. 
16); the distribution extends down to Ltir/ Lr ^ 1 with a 
sudden drop of the galaxy number at L^m /Lr < 1, indicat- 
ing a paucity of optically thin galaxies (t^"*^ ^1)- This is 
also indicated by the sequence of Mdm/Mz{iii) oc ^^(^hi) 
in Fig. 25. If there are a significant number of opti- 
cally thin galaxies, such galaxies should appear in the 
region of ^^(hi) ^ 1 Mq pc~^ with a constant value 
of Mdm/Mz(Hi) ^ 1, corresponding to Ed < E^^crit and 
Mdin/Md ^ 1- However, such galaxies are not found, and 
the sequence seems to stop at ^ E^^crit {t^°cS 1)- 

If the optical thickness to dust-heating radiation does 
not have any effect on galaxy formation efficiency, a fine- 
tuning would be required for the sharp drop of galaxy 
distribution around T^°^ff ^1. A possible explanation 
is a selection effect; optically thin galaxies would have 
small infrared-to-optical luminosity ratio and hence could 
be missed in the samples selected by infrared emission. 
To check this possibility, we plot in Fig. 27 velocity- 
integrated 21-cm flux density (5int) versus angular size 
(Opgo) for the SDSS-Hl and HIPASS samples, in compar- 
ison with the AKARI-SDSS and AKARI-HIPASS galax- 
ies. The lines corresponding to several values of H I col- 
umn density Eri = Mh_j/{ttR'^) are also depicted. It can 
be seen that there is no significant difference between 
the distributions of Eri for the AKARI-detected and - 
nondetected galaxies, and the paucity of galaxies in the 
region of Ehi ^ 10 M© pc~^ can be seen, roughly corre- 
sponding to T^°*ff ^ 1 for typical values of /hi, fdz, and 
Z. Hence the infrared selection bias cannot explain the 
paucity of optically thin galaxies. It should be noted that 
the beam size of HIPASS is 15.5 arcmin, which is much 
larger than the galaxy sizes, and hence there should be no 
selection bias about surface brightness in the H I data. 

Therefore we consider that there is a physical effect re- 
ducing the efficiency of galaxy formation when rjf^g ^ 1. 
Because the opacity is to dust-heating radiation, it is most 
likely the radiative feedback by the dust photoelectric 
heating of ISM. A part (typically epE 1-10%) of ra- 
diation energy absorbed by dust grains is converted to 
photoelectrons and subsequently heat ISM, and this ef- 
fect is large enough to be the dominant heating process 
in many phases of ISM in galaxies (see, e.g., Bakes & 




Fig. 26. The properties of galaxies in the AKARI-SB sample. The upper-left, upper-right, and lower-right panels are the same as 
Figs. 18, 21, and 22, respectively, but showing the AKARI-SB galaxies by the filled red squares, in comparison with the AKARI-SDSS 
sample (open triangles). The lower-left panel shows the histogram of the mass ratio of heated dust to hydrogen gas, Mjjjr/Mjj, 
where Mh is H2 and Hi mass for the AKARI-SB and AKARI-SDSS samples, respectively. The histogram of the AKARI-SDSS 
sample is multiplied by a factor of 0.1 for comparison. 



Tielens 1994; Wolfire et al. 1995; Rubin et al. 2009). 
Then, in T^°^g ^ 1 galaxies, it is expected that the lack of 
self-shielding of dust-heating radiation leads to the sup- 
pression of star formation activity and reduction of neutral 
gas in ISM on a galactic scale. Bigiel et al. (2010) found 
that star formation efficiency becomes extremely low with 
decreasing Shi in outer disk of nearby galaxies, which is 
consistent with the picture discussed here. 

On the other hand, in the large limit of Ehi, Fig. 27 
shows the paucity of galaxies at the region of Ehi ^ 
10^ Mqpc"^. This is likely because most of gas tends 
to be molecular when gas column density is high, result- 
ing in a saturation of Ehi, as observed for nearby galaxies 
(Bigiel et al. 2008). The efficient self-shielding of dust- 
heating radiation may be responsible for converting H I 
into H2. The lack of both large and small Ehi galaxies by 
these effects may explain the known good correlation be- 
tween Hi mass and size, Mhi oc (Haynes & Giovanelli 
1984; Rosenberg & Schneider 2003; Garcia-Appadoo et al. 
2009; van der Kruit & Freeman 2011). The dispersion in 



Mhi/R^ seems especially small when Hi size is used for 
R (Rosenberg & Schneider 2003), further indicating that 
rjf^ff is playing an essential role. 

6.2. Implications for the Kennicutt- Schmidt Relation 

The quantities about infrared emission discussed in this 
work are mostly in surface densities of galaxies, e.g., ?7tir 
or SdiR- These quantities are closely related to the sur- 
face density of SFR, Ssfr, and surface hydrogen gas den- 
sity, Sh- Since the correlation between Esfr and Eh is 
frequently discussed in the literature as the Kennicutt- 
Schmidt relation (Schmidt 1959; Kennicutt 1998), it is 
interesting to discuss about it in light of the results ob- 
tained in this work. 

In Figure 28 we show the Esfr-Eh plot for the AKARI- 
SDSS and SDSS-Hl samples, using Hi gas mass for Eh- 
For the K98 SB sample, we plot the original data for the 
36 galaxies in K98 (Eh from H2 mass by CO observa- 
tions), corrected only for SFR into the Kroupa IMF. The 
SFRs of the SDSS galaxies are from optical spectra, while 
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angular size dpgg [arcsec] angular size dpgg [arcsec] 

Fig. 27. Velocity-integrated flux density of the 21-cm emission (Sint) versus angular size (Spgo), for the AKARl-SDSS (left) and 
AKARI-HIPASS (right) samples. The grey dots are for galaxies whose infrared emission was not detected by AKARI. The solid 
lines correspond to several values of Hi column density Shi as indicated in the panels (Sint oc Mm/d^, 9pqq = R/dA, and hence 
Eh I oc Sint /'^pgo ) ■ For these solid lines we assumed a redshift z = 0.02 that is typical for the samples here, though the dependence 
on z is very small only by the cosmological effect proportional to (1 + z). 



those of the K98 SB sample are from Xtir- There is no 
good correlation between Ssfr and Shi for galaxies in 
the SDSS-based samples at Eh ^ 100 Mopc"^, which is 
consistent with the results of K98 for normal spirals and of 
Bigiel et al. (2008) for sub-kpc regions of nearby galaxies. 
This range of Eh roughly corresponds to t^^^ ^ 1, and as 
argued in §6.1, the radiative feedback process seems to be 
working to suppress the star formation activity. Then the 
relation between SFR and Hi gas amount would become 
complicated by unstable competition between star forma- 
tion and its feedback. This may give an explanation for 
the large dispersion in the Esfr-Ehi relation. 

On the other hand, we have found in §5.5 that Uh is 
constant at Uh ~ 50 for the AKARI-SB sample, corre- 
sponding to a constant SFE of i/'e,max ^ 10 Gyr^^, i.e., 
^SFR cx Eh. This SFE is equivalent to the star formation 
time scale of isF,min = V'tTmax 0.1 Gyr. It is interesting 
to note that this characteristic SFE is in nice agreement 
with those found in the correlation between Mh2 and SFR 
for a wide variety of galaxies as well as molecular clouds 
in the Galaxy (Gao & Solomon 2004; Wu et al. 2005; 
Bigiel et al. 2008; Evans et al. 2009), when the gas mass 
is measured from tracers of dense molecular gas. These 
results indicate that i/'e.max is a universal characteristic 
value originating from the physics of star formation in 
ISM, rather than time scales on a galactic scale. It is 
widely known that the efHciency of star formation from 
molecular clouds is typically a few percent in mass frac- 
tion, and typical dynamical time scales of molecular clouds 
are Myr, resulting in a star formation time scale of ~ 0.1 
Gyr. This maximum SFE ■0e,max seems to be achieved on 
a galactic scale when t^^^^ > 1 (Eh > 100 Mq pc~^) be- 
cause the self-shielding effect of dust-heating radiation is 
working for efficient gas cooling. Most of interstellar gas 
would be molecular in such galaxies. 

Now we may interpret the overall trend of the galaxy 



distribution on the EgpR-En plane by the strong depen- 
dence of SFE on the dust opacity, as follows. There 
are very few galaxies having t^Iq ^ 1, because the feed- 
back suppresses an efficient star formation. Galaxies ap- 
pear in the region of r'^^ff ^ 1 with a large dispersion 
in the Kennicutt-Schmidt law, because SFE rapidly in- 
creases around T^°*ff ^ 1. Then SFE continues to in- 
crease and asymptotically reaches the maximum value of 
V'e,max ^ 10 Gyr~^ at T^°eff » 1. We expect many galax- 
ies around T^°*ff ^ 1 because star formation time scale is 
relatively longer than that in more dusty galaxies. 

6.3. Implications for Theoretical Modelings of Galaxy 
Formation 

There are a few implications of the results obtained by 
this work for theoretical modeling of galaxy formation and 
evolution. First of all, the tight correlation between Uh 
and J7tir for normal galaxies provides a convenient way to 
calculate infrared SEDs in cosmological galaxy formation 
models, in which the total dust-unabsorbcd luminosity of 
stellar emission and size of a galaxy can be calculated in 
a relatively straightforward manner. When galaxies are 
optically thick in terms of T^°gfj as indicated by our study, 
the total stellar luminosity gives a direct estimate of Ltir. 
Then Uh or can easily be calculated from the C/tir- 
Uh relation, which is more reliable than LTiR-Td relation 
because of the tighter correlation and more solid physical 
background. There is no particular reason to expect that 
this relation does not hold for high redshift galaxies, if 
the dust properties are not significantly different and the 
spatial distribution of stars and dust grains is the case 
B2. It must be examined whether the existing ah initio 
galaxy formation models of infrared emission satisfy this 
relation, especially at z ~ 0. 

Another important implication is the strong depen- 




Fig. 28. SFR column density versus hydrogen gas column 
density (the Kennicutt-Schmidt law). The red squares are for 
the circumnuclear starbursts of the K98 sample, estimating 
gas mass from H2 and SFR from Ltir- The AKARI-SDSS 
galaxies are plotted by triangles, and grey dots are for SDSS 
galaxies without AKARI detection, estimating gas mass from 
H I and SFR from SDSS spectra. The relation corresponding 
to star formation efficiency (SFE) = 10 Gyr^-' is shown by 
the solid line. 

dence of star forming activity on dust opacity or gas col- 
umn density suggested in §6.1 and §6.2. Various feedback 
processes have intensively been discussed in the context of 
cosmological galaxy formation (e.g., Benson et al. 2003; 
Baugh 2006; Silk 2011), but these feedback processes are 
not physically connected to the galaxy-scale dust opacity, 
T^°^*jj. Photoionization of hydrogens may have a similar 
dependence on gas column density, but it has a much 
larger opacity than r^^ff by a factor of ^ 3 x 10"* for 
Md/Mn = 0.01. Supcrnova/AGN feedbacks are unlikely 
related to the galaxy-scale dust opacity. 

We have argued in §6.1 that the dust photoelectric heat- 
ing of ISM is the likely process responsible for the feed- 
back working in T^°^ff < 1 galaxies. Using eq. (4), 1 M© of 
star formation produces 3.2 x 10""'(epE/0.03) erg by photo- 
electric heating, when most of heating radiation energy is 
absorbed into dust grains. This heating energy is compa- 
rable with and may even be larger than 1.1 x lO**^ erg pro- 
duced by the supernova feedback (assuming 10^^ erg per 
one supernova, 8 M© threshold for supernova explosions, 
and the Kroupa IMF). In contrast to supernova feedback, 
which is a local energy input only to dynamically con- 
nected regions, photoelectric heating could suppress star 
formation on the global scale in a galaxy when t^"^^ ^ 1. 

Such a feedback effect as a function of dust opacity 
(or gas column density) may systematically change the 
galaxy formation and evolution. Cosmological objects 
with dark matter mass Mdm virializing at redshift z have 
virial density proportional to pvir oc (1 -l-z)'^, and the avail- 
able baryon mass and galaxy size roughly scale with dark 
matter mass and virial radius, respectively. Therefore we 
expect that gas column density scales as oc M^^^{1 + z)^. 
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indicating that this effect will accelerate star formation 
at higher redshifts, also weakly favoring star formation in 
massive objects at a fixed redshift. 

This may have interesting implications for some ma- 
jor problems in the hierarchical galaxy formation scenario 
(e.g., Baugh 2006; Silk 2011). One is the difficulty to 
reproduce the early appearance of massive and quiescent 
galaxies at high redshifts (Cimatti et al. 2004; Glazebrook 
et al. 2004), and another is that the theory predicts too 
large stellar mass in massive dark halos 10 ^^M©) at 
z ^ compared with the observed galaxy luminosity func- 
tion (the so-called "ovcrcooling" problem). If SFE is a 
strong function of T^°^ff, it would increase the star forma- 
tion activity at high-redshift compact dark halos. Massive 
objects have larger gas column density at a fixed redshift, 
which may be relevant to reproduce the well-known trend 
of downsizing in galaxy formation. On the other hand, 
dark halos collapsing at low redshifts are massive but have 
low densities, and star formation would be suppressed to 
alleviate the overcooling problem. 

6.4- Comparison with Previous Studies 

One of the most important results of this paper is 
the tight correlation between the dust-heating radia- 
tion strength Uh and the galaxy-scale infrared radiation 
strength C/tir- The parameter [/tir is equivalent to the 
infrared surface brightness, and the correlation of this 
quantity to infrared SED has been studied by Lehnert & 
Heckman (1996, hereafter LH96) and Chanial et al. (2007, 
hereafter C07). 

LH96 estimated the size of 32 nearby starburst galax- 
ies by Ha images, and found a correlation between 
Ltir/ [T^Rna)'^ and infrared color temperature estimated 
by 60/100 fim flux ratios. They argued that a surround- 
ing dust screen radiated by a radiation field strength of 
-^tir/(7''^Hq)^ is an adequate zeroth-order description of 
their data, which is equivalent to our interpretation of the 
case B2 for normal star-forming galaxies. However, we 
also found that Uh becomes roughly constant for circum- 
nuclear starbursts in the high SFR density region. The 
LH96 galaxies are in the range of 10^ ^ Ltiji/{ttR'^^) ^ 
10^ Lq pc~^, roughly corresponding to 10 ^ Utir ^ 10'^ 
assuming that the Ha radius is not different from the r- 
band radius used here. This is the region of transition 
from Uh oc ?7tir to a constant Uh, and indeed, the distri- 
bution of the LH96 galaxies is fiattcr than expected from 
the simple surrounding screen in the Seofim/ SiQOfj,^ ver- 
sus LTiR/(7ri?^^) plane (see their Fig. 5). Our result is 
thus not inconsistent with that of LH96, but revealed the 
correlation more quantitatively by larger statistics, using 
not only starbursts but also normal galaxies. 

C07 studied on samples (~400 galaxies in total) in- 
cluding non-starburst galaxies. They estimated galaxy 
sizes by radio continuum (RC) images, and found a cor- 
relation of Stir = Ltir/{2t:RI^) oc with S - 18.5, 
where is the dust temperature estimated by assum- 
ing the modified blackbody spectrum. The range of Stir 
is ~ lO^-lO^L© pc-2 {Utir ^ 0.1-10^). In contrast to 
LH96, they concluded that the surrounding dust screen 
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is ruled out because the observed slope of the correlation 
is significantly different from the expectation ((5 = 4 + /3). 
Instead, they proposed another interpretation based on 
the Kcnnicutt- Schmidt law. Assuming a constant dust-to- 
gas ratio CSdiR oc Sgas), -^tir oc -0, and the original form 
of the Kennicutt- Schmidt law (Ssfr ^gati Kennicutt 
1998), the basic equation of Ltir oc M^irT^"''^ can be re- 
duced to Stir k y3.5(4+,9)^ which is close to the observed 
6 value with the emissivity index of /3 ~ 1-2. 

Our result of ?7tir <x Uh (x T^^^ for normal galaxies 
(t'^TiR ^ 30) is equivalent to 5 = 4-1-/3, while the constant 
Uh for the circumnuclcar starbursts ([/tir ^ 30) corre- 
sponds to 5 = 00. The 6 value obtained by C07 may be a 
result of the fit including the transition region of J/tir be- 
tween the two regimes. Their interpretation based on the 
Kennicutt-Schmidt law is not supported from our study, 
because wc have seen that SFR is not proportional to 
Xtir, and S 

diR is roughly constant against change of Sgas, 
for normal star-forming galaxies. There is a large scatter 
in the Kennicutt-Schmidt law, and it is unlikely that the 
tight correlation between Uh and ?7tir emerges from this 
law. 

Besides the larger statistics, the strength of our study 
is that we estimated the dust-heating radiation field 
strength by fitting realistic physical dust models to the 
four AKARI/FIS photometric bands at 65-160 cov- 
ering the broad thermal SED peak. It is good to use Uh 
instead of Td to examine the correlation against C/tir, be- 
cause they have the same physical dimensions and hence 
not only the scaling between the two but also the ratio 
Uh/Utik is physically meaningful. It should also be noted 
that Td is meaningful only for large grains in thermal equi- 
librium. Another important source of difference is that, 
as mentioned in §5.2, the AKARI-bascd samples include 
galaxies with smaller Uh than the IRAS-based samples 
used in the previous studies, which enabled us to clearly 
find the linear correlation between Uh and C/tir in the 
normal galaxy regime. 

7. Summary and Conclusions 

We have constructed two new low-redshift galaxy sam- 
ples based on the AKARI/FIS All Sky Bright Source 
Catalog: the AKARI-SDSS sample (878 galaxies) cross- 
correlated with the SDSS galaxies and available Hi data, 
and the AKARI-HIPASS sample (711 galaxies) cross- 
correlated with the HIPASS H I galaxy catalog. To com- 
plement these for wider dynamic range of surface SFR 
density, the AKARI-SB sample (24 galaxies) has also 
been constructed from nearby circumnuclcar starbursts. 
These samples include only galaxies having high quality 
AKARI/FIS fiuxes at least in two bands to ensure the 
reliability of SED fittings. 

We then studied these samples to understand the 
physics determining galactic-scale infrared SEDs, espe- 
cially the broad thermal peak by large grains. The 
strength of dust-heating radiation field {Uh, roughly corre- 
sponding to the dust temperature as Uh oc T'| ''^), heated 



dust mass (Mdm), and total infrared luminosity (Ltir) 
were derived from the fittings of physical dust models to 
the observed infrared SEDs, and their correlations with 
various physical quantities [optical size and luminosity, 
SFR {ip), stellar mass (Mh.), metallicity, and Hi or H2 gas 
masses] were investigated. 

We found a tight power-law correlation between the 
^tir/V' ratio and the specific SFR (SSFR, ips = il^/A'h), 
as LTm/4' oc "07 "'^^j which is valid except for very dusty 
galaxies having a large infrared-to-optical luminosity ratio 
(§5.1). Therefore SSFR is a good measure of the contribu- 
tion from young stars to Ltir- This relation may be use- 
ful to estimate the third parameter when two of SFR, M^, , 
and £tir are known. The widely used value of Ltib./iP as 
a SFR indicator is valid only for galaxies having largest 
SSFR of -03 ^ 0.1 Gyr~^, and relatively aged stellar pop- 
ulation must be taken into account in the dust-heating 
energy budget for most of the AKARI-SDSS galaxies. 

The most important result of this work is the discovery 
of the tight linear correlation between Uh and galactic- 
scale mean infrared radiation field C/tir oc £tir/ {t^R^) for 
the AKARI-SDSS and AKARI-HIPASS samples, taking 
the r-band optical size as R (§5.3). The ratio of C//i/C/tir 
is close to unity in the same physical units, and the phys- 
ical dispersion along this relation in Uh is about 0.3 dex 
(^ 13% in the dust temperature). This C/tir-CZ/i relation 
is nicely explained by a thin disk of dust-heating sources 
embedded in a thicker, optically-thick {t^^q ^ 1) dust 
disk (the case B2 in Fig. 15), where T^lg is the effective 
dust opacity averaged over wavelength with the weight of 
the SED of the interstellar radiation field (ISRF). This 
also means that infrared SEDs of galaxies are determined 
mainly by ISRF on a galactic scale, rather than those 
on smaller scales like individual star forming regions or 
molecular clouds. 

The C/tir-CZ/i relation is particularly useful in theoreti- 
cal modeling of infrared emission, giving a simple method 
to determine dust SED from Ltir and R that can be pre- 
dicted relatively easily in cosmological galaxy formation 
models. It may also be useful in observational studies; 
e.g., when flux of dust emission is available only in one 
band for a galaxy but its size is known, one can make a 
guess about the SED shape and Ltir by searching for Uh 
that satisfies the relation of C/tir ~ Uh- 

Since only dust in the layer of one optical depth is 
heated, the dust column density estimated by infrared 
emission should be constant at S^ir ~ S^.crit correspond- 
ing to Td^oS = 1, which was also confirmed by our data. 
This interpretation predicts that the heated fraction of 
the total dust mass {Mdui/^Id) scales with the total dust 
column density as cx S^^, and the observed trends of 
MdiR/M2(Hi) cx S^^^jjjj and A/diR/Mni oc Sjjj support this 
expectation, where 2'(Hl) denotes the metal in Hi gas 
phase (§5.4). An important implication is that, for many 
normal galaxies, the dust mass derived from infrared emis- 
sion gives only a lower bound, with a significant amount 
of cold dust not contributing to infrared emission. 

On the other hand, the data of the AKARI-SB sam- 
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pie indicate that there is an upper limit of [//i ^ 50 (in 
units of the local ISRF around the solar neighbourhood) , 
and Uh becomes constant at this value for intensive star- 
bursts having Utir ^ 50 (§5.5). The H2 mass is tightly 
correlated with the heated dust mass M^ir, and the upper 
limit of Ufi corresponds to the maximum star formation 
efficiency (SFE, defined as SFR per hydrogen gas mass) of 
■06, max ^ 10 Gyr~^. This maximum SFE likely originates 
from the physics of star formation in ISM, and the spatial 
distribution of heating sources and dust changes from the 
case B2 to Bl (the same mixed distribution for stars and 
dust) in Fig. 15 around J7tir ^ 50. 

All the galaxies studied here seem optically thick to 
dust heating radiation. The distribution of T^°^g extends 
down to T^°eff ~ ^' there is a sharp drop of the number 
of optically thin galaxies at T^°^ff < 1. We argued in §6.1 
that this is unlikely to be a result of selection effects, but 
there seems a feedback effect to suppress the formation 
of optically thin galaxies. Because it is related to dust 
opacity, the feedback is likely by the photoelectric heat- 
ing of ISM by dust grains. This process is known to be 
the major heating in many phases of ISM, and the self- 
shielding of dust-heating radiation on a galactic scale may 
be a necessary condition to avoid a global suppression of 
star formation. It is often stated that starburst galaxies 
are dusty, but if this hypothesis is true, one may rather 
state that they are starbursts because they are dusty. 

We discussed implications of our results for galaxy for- 
mation in general. The strong dependence of galactic- 
scale gas cooling on T^°*g may be responsible for the known 
correlation between H I mass and size (Mhi (xR^), which 
is on the line corresponding to T^°gg ~ 1 (§6.1). The 
Kennicutt-Schmidt law (SFR surface density Ssfr versus 
hydrogen gas surface density Eh) may be understood by 
the maximum SFE V'e.max as the baseline with the strong 
dependence of SFE on T*°*fj by the feedback (§6.2). The 
feedback by dust opacity may also have interesting impli- 
cations for galaxy formation in the cosmological context, 
because it would accelerate star formation in massive ob- 
jects at high redshifts. Such an effect may be helpful to 
solve some of the problems currently discussed in the field 
of cosmological galaxy formation (§6.3). 

It is quite intriguing to examine whether the scaling 
relations found in this work, especially UTiR-Uh, hold also 
for high redshift galaxies. Planned infrared/submillimeter 
facilities will give such opportunities in the near future. 
Especially, ALMA will be a quite powerful facility for this 
purpose by the unprecedented flux sensitivity and angular 
resolution at wavelengths covering the thermal SED peak 
at high redshifts. 

The electronic catalog data of the samples constructed 
in this work are available on request to the authors. 
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